The Double Exchange Model at Low Densities 
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We obtain the phase diagram of the double- exchange model at low electronic densities in the pres- 
ence of electron-electron interactions. The single particle problem and its extension to low electronic 
densities, when a Wigner crystal of magnetic polarons is generated due to unscreened Coulomb inter- 
actions, is studied. It is argued that the Wigner crystal is the natural alternative to phase separation 
when the Coulomb interaction is taken into account. We discuss the thermal and quantum stability 
of the crystalline phase towards a polaronic Fermi liquid and a homogeneous, metallic, ferromagnetic 
phase. The relevance and application of our results to EuBe is also considered. 
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PACS numbers: Pacs: 71.10.-w,73.20.Mf,73.20.Qt,75.47.Gk 

Introduction. Magnetic polarons are ubiquitous in the 
physics of double-exchange systems such as mixed valence 
manganites 0, antiferromagnetic and ferromagnetic 
semiconductors (such as EuO 3]), diluted magnetic semi- 
conductors (DMS) (such as Gai-^Mn^As) 0,01, and in 
EuB 6 where they have been observed in the optical re- 
sponse 0. Colossal magnetoresistance (CMR) has been 
attributed in part to the presence of magnetic polarons 
close to the paramagnetic- ferromagnetic transition [7| . In 
fact, CMR in the Mn pyrochlores has been proposed on 
the basis of magnetic polarons [j| . Given that these dif- 
ferent classes of magnetic materials have attracted con- 
siderable attention in recent years because of their po- 
tential for the development of new magnetoelectronic de- 
vices, a number of theoretical approaches to the problem 
have been developed through the years. In particular, ex- 



tensive work has been done with em pha sis in the physics 
of magnetic semiconductors iilUa and CMR man- 
ganites [ri| . where electron-electron interactions are as- 
sumed to be weak. 

In this work we will focus our attention on the phase di- 
agram of the double-exchange model (DEM) at low den- 
sities, taking into account the effects of Coulomb interac- 
tions between the electrons. Our results are summarized 
in Fig. ^ Studies devoted to the polaronic stability in 
this particular model and its variations have been per- 
formed by several authors, both analytically and numeri- 
cally 01 ) but no unified theory has been presented so far. 
Here we study the extreme case of zero density (that is, 
the single magnetic polaron), the Wigner crystallization 
of magnetic polarons due to Coulomb interactions, and 
its thermal and quantum melting into a polaronic Fermi 
liquid and a ferromagnetic metal. We argue that the 
phase separation instability of the non-interacting DEM 
at low densities is replaced by Wigner crystallization and 
a polaronic Fermi liquid, stabilized by electrostatic in- 
teractions. We also show that the mean-field approach 
breaks down in the low density regime and that electron- 
electron interactions have to be taken into account prop- 
erly. 

The Model. In the presence of the Coulomb term, the 
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Figure 1: Phase diagram of the diluted DEM model as a 
function of temperature and density (t — 0.5 eV, a = 4 A): 
paramagnetic (PM) , polaronic (Pol) , polaronic Wigner crystal 
(WS), and ferromagnetic (FM). 



DEM reads (we use units such that h = 1 = ks)' 
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where ai (<zj) is the annihilation (creation) opera- 
tor for an electron on site r; with its spin aligned 
in the direction of a classical localized spin Si = 
S(sin0icos</>j,sin0isin<^j,cos0j), Oi and <j>i being the lo- 
cal Euler angles. In Eq. QJ, t tj — t[cos {9 l /2) cos (6j/2) + 
sm(0j/2)sm(0j/2)e -l ^ i ~'^] is the hopping energy, e 
the electric charge, rij = a\a,i the number operator, and 
n e the average number of electrons per atom. This model 
can be obtained from an interacting Kondo lattice Hamil- 
tonian in the limit where the exchange interaction be- 
tween local and electron spin, Jh, satisfies Jh ^> t 0. 

The presence of the full Coulomb interaction in eq. (JTJ 
is to be understood in connection with the case of ex- 
tremely reduced electron density. In a conventional Fermi 



2 



0.4 




Figure 2: Top: DOS for the 2-D spin configuration shown on 
the bottom (exact diagonalization, bars) and averaged DOS 
for the homogeneous case (recursive method, continuous line). 
The magnetic polaron level (marked by an arrow) agrees with 
— 4icos[7r/(ii + 1)]. Bottom: Magnetic polaron wavefunction 
(top) and its underlying lattice spin configuration (bottom). 



liquid the high density of electrons makes the screen- 
ing process very effective, and the effect of the electron- 
electron interactions can be absorbed into the renormal- 
ization of physical quantities such as the effective mass. 
In a very diluted electron gas the Coulomb interaction 
cannot be addressed meaningfully in this way. In the 
language of the one component plasma, this can be un- 
derstood with reference to the dimensionless parameter 
T s = r /ao, where r Q — (S/Aitp) 1 ^ (p being the volumet- 
ric density) is the average distance between electrons, and 
ao is the Bohr radius. While the kinetic energy scales as 
l/tg, the potential energy varies as l/r s , and, therefore 
dominates in the low density (r s — > oo) regime. It is well 
known since Wigner |l3j | that, under such circumstances, 
the electrons arrange themselves in a regular lattice. An- 
other aspect related to the reduced carrier density is pre- 
cisely the stability of free magnetic polarons. Within the 
simple DEM, the only possibility for the stabilization of 
magnetic polarons is at low electronic densities because, 
the lattice spins of several neighboring unit cells are ex- 
pected to participate in the magnetization cloud of each 
electron. Were it otherwise, the electronic wave functions 
would overlap considerably, destroying the polaron pic- 
ture. On the other hand, if a Wigner crystal is realizable 
in a system described by eq. then, as in any solid, 
there is zero point motion of the electrons about their 
equilibrium positions. If the electron itinerates among 
several unit cells during this motion, the magnetic cou- 



pling, Jh, leads to the local polarization of the lattice 
spins, producing a bound magnetic polaron. This pro- 
cess is obviously limited by the melting of the electronic 
solid and, thus, the question arises of how to describe 
these two tendencies for the polaron formation and the 
interplay of Coulomb and magnetic interactions. A final 
aspect that cannot be neglected in this regime is the fact 
that the classical approximation for the spins implies an 
intrinsically disordered electronic problem, at any non- 
zero temperature (a simple consequence of the depen- 
dence of tij upon the local spin texture). This requires, 
for instance, the calculation of the spectral properties for 
a disordered electronic s yste m, executed here by means 
of the recursion method [lj| and exact diagonalization. 

Single Particle Regime. At zero thermodynamic den- 
sity the Coulomb interaction can be disregarded, and the 
problem reduces to a single magnetic polaron. We con- 
sider the electronic density of states (DOS) for a given 
spin configuration, {S^}, generated with a Boltzmann 
weight that includes the mean-field required to obtain 
a given average magnetization, M, of the lattice spins 
|1(t| |. Averaging over several realizations of disorder, we 
can write the electronic energy as: 

E el (M,n e ) = J G(E F {n e ,M)-e)e{N(e,M))de, (2) 

where the explicit dependence of both the DOS (N(e)) 
and Fermi energy [E F ) on the magnetization is empha- 
sized. For n e <C 1 we approximate eq. (J2J simply by 
E e i(M,n e ) ss Eb(M)n e , with Et representing the bot- 
tom of the band that, in 3-D, is found to lie at — At for 
the paramagnet (Ai — 0), and at — 6t in the fully po- 
larized ferromagnetic phase (Ai = 1). Let us assume 
that the polaron consists of a cubic volume of side R (in 
units of the lattice parameter, a), inside which Ai = 1. 
Outside the polaron the system remains paramagnetic 
(Ai =0). In Fig. [2] we show exact diagonalization re- 
sults of a two-dimensional ferromagnetic region embed- 
ded in a paramagnetic phase, leading to the formation 
of a bound state, that is, a magnetic polaron 17]. The 
free energy per electron, at a temperature T, required to 
form a magnetic polaron from the paramagnetic phase is 
given by: 

AF Pol (R, T) /n e = 4t~6tcos(Tr/(R+l)) 

+TR 3 \og(2S + l)-TS C f g (n e ,R), (3) 

where 5cfg is the configurational entropy. The first three 
terms in reflect the two competing effects at play: the 
first is the electron's preference for a ferromagnetic back- 
ground, accompanied by an energy cost for localization; 
the second is the reduction of entropy caused by the ap- 
pearance of the (fully polarized) magnetic polarons. For 
n e <C 1 we find that Scf g (n e , R) is always negligible (l7| . 
Minimizing eq. J3J {R ^> 1) with respect to R gives: 

R cq (T) ~ [2t7T 2 /(Tlog (2S + 1))] 1/5 (4) 
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for the equilibrium radius of the magnetic polaron, in- 
creasing at low temperatures as T -1 / 5 . The stability 
condition for the polaron is AFp Q \(R eq ,T) < 0, which 
provides the temperature, T m , below which the polarons 
first appear: 

T m /t = 8\/2/(25v / 57r 3 log(25 + 1)) . (5) 

Note that at T = T m the polarons are formed with a 
size given by R C q{T m ), and hence, the crossover from the 
paramagnetic to the polaronic phase is discontinuous. 

The previous results can also provide an estimate of the 
temperature below which a homogeneous ferromagnetic 
state is established. Assume that each electron forms 
a polaron that is randomly distributed over the system. 
Since the polaron radius increases according to as the 
temperature is lowered, there is a critical temperature, 
Tq , below which the polarons percolate into a ferromag- 
netic phase. This percolation threshold is defined from 
the criterion n e R^ q ~ p c (p c sa 0.31 for a cubic lattice), 
yielding from Q a Curie temperature, 

Tg/t =i 27r 2 /[log (25 + 1)] (n e /p c f 3 . (6) 

This temperature can be compared with the Curie tem- 
perature obtained by a mean- field approach : T c SF /t = 
5(1 + S)/(3S)n e . Notice that, for n e — ► 0, we have 
Tq <C Tq IF , indicating that the mean-field calculation 
overestimates the Curie temperature in the low density 
regime. For Tc < T < T m , an anomaly in the paramag- 
netic susceptibility is expected to signal the polaron pres- 
ence, through an enhanced effective moment. Naively, 
the transition from the polaronic phase to the ferromag- 
netic phase would be first order, for the magnetization 
jumps at Tc- Nevertheless, because this transition is 
percolative in nature, one would expect a second order 
phase transition, in the sense that the average magne- 
tization should be weighted by the "mass of the infinite 
cluster", which evolves continuously from the percolation 
threshold. 

Polaronic Wigner Crystal. The previous estimate of 
the Curie temperature for a small but finite density of 
electrons has a fundamental problem: the non-interacting 
DEM is unstable towar ds p hase separation at low densi- 
ties E E3 EJ |H . In the presence of Coulomb 
interactions, however, the phase separation is frustrated 
by the large electrostatic energy price of confining charge 
to a given region of the system. It is found that, for 
reasonable values of the dielectric constant, the tendency 
for charge neutrality, which favors small radius of elec- 
tron rich regions, and the kinetic energy of localization, 
lead to a strong suppression of the phase separation re- 
gion, in temperature and electron concentration, and its 
replacement by a polaronic phase [l7j . 

At the lowest densities (r s 3> 1) the Coulomb inter- 
action is unscreened and larger than the kinetic energy 
in 0. In order to minimize the Coulomb energy, the 
electrons form a Wigner crystal and undergo zero point 
motion around their equilibrium position. Derealization 



of the electron wavefunction due to quantum fluctuations 
leads to the polarization of the local spins around the 
crystalline positions. The local polarization, exactly as 
in the case of the single polaron discussed previously, pro- 
duces an extra confining potential to the electron. Fol- 
lowing the Wigner-Seitz approximation |15| . the Wigner 
crystal unit cell (much larger than the original lattice 
spacing a) is approximated by an electrically neutral 
spherical volume, inside which the ionic charge density 
is homogeneous. The electrostatic potential energy then 
depends only upon r: the distance of the electron from 
the center of the cell. The Hamiltonian for an electron 
in this uniform charge and spin background is then: 

where p is the electron momentum, m = l/(2a 2 t) the 
effective electron mass, and uj 2 = w|/3 = e 2 /mr 3 , where 
uip is the plasma frequency. In Eq.(7jl u is the frequency 
of the confining potential due to the DEM mechanism, 
and is a variational parameter. The radius of the mag- 
netic polaron in the Wigner crystal relates to lu by: 

R = ^itjn, (8) 

where fl = \Ju) 2 + u> 2 is the total frequency of oscillation 
of the electron. Notice that the ground state energy of 
is E a = -6t-3e 2 /(2r ) + 3fi/2 and hence the relative 
gain in free energy is: 

ATwp = -2t + 3 (fi - uo e ) /2 + 4irR 3 T log (25 + 1) /3. 

(9) 

Minimization of eq. © with respect to u> gives: 

R(T) = { gCrV*) 1 ' 5 , (io) 

where Rs — [3t/(w e )] 1 / 2 is the saturation radius, and 
T*/t = 9/[47i\R|log(25+ 1)] is the temperature below 
which the polaron radius saturates due to the interplay 
between the DEM and the Coulomb interaction. 

Wigner Crystal Melting. It is clear that the previous 
results are valid for temperatures so low as not to melt 
the Wigner solid. The calculations for the independent 
polaron model reveal that the temperatures for polaron 
stability are already typically small for reasonable val- 
ues of t, but the electronic solid is much more sensitive 
to the temperature. The Wigner crystal melting tem- 
perature, Tm, can be estimated from the Lindemann's 
criteria j2(J: Tm w 0.01(e 2 /a)rh/ 3 . It is known from sev- 
eral numerical calculations HM Eil on the stability of 
the one component plasma, that the maximum densities 
and temperatures at which the Wigner crystal can exist 
correspond to r s ~ 50-100, and T ~ 10 K. Values of 
r s ~ 50-100 correspond to n e - lO^-lO" 5 for t = leV 
and a = 4 A. Because of the absence of magnetic inter- 
actions between different polarons, the Wigner crystal is 
a superparamagnet. In the presence of other long-range 
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interactions (such as dipole-dipole) the polaronic Wigner 
crystal can exibith long range magnetic order. Increasing 
the electron density at T = causes the Wigner crystal 
to quantum melt at some critical density with two possi- 
ble outcomes: a paramagnetic polaronic Fermi liquid or a 
fully polarized ferromagnet. In both cases the carriers are 
mobile and can screen the long-range part of the Coulomb 
interaction leading to a Fermi liquid state. At finite tem- 
peratures, where the electron state cannot be described 
by the zero point motions implicit in eq. 10 alone, the 
crystal should follow the features of the phase diagram 
for the electron gas j2(| . The characterization of the sys- 
tem in the neighborhood of the melting point, where the 
presence of a polaron liquid is plausible (Fig.Q), is restric- 
tively hard, even for the simple electron gas. Far from 
this region, where the electron density is high enough to 
make the screening process effective, one expects to re- 
trieve the behavior obtained before within the indepen- 
dent polaron model discussed previously. Our results are 
summarized in the phase diagram presented in Fig. ^ 

Application to EuBq. While we have produced the 
general phase diagram for the DEM model at low den- 
sities, it is interesting to apply our results to a real sys- 
tem. EuBg is a good magnetic metal with extremely 
reduced electron density [221, exhibiting all the charac- 
teristic signatures of a polaronic phase in Raman scatter- 
ing . Such experiments reveal that the FM transition 
at T c * p ~ 15 K is preceded by an interval of temper- 
atures (T < T°f v as 30.fr) where the system is domi- 
nated by the presence of magnetic polarons. We have re- 
cently proposed that EuBg is a DEM material in the low 
density regime [l6[ . characterized by a hopping integral 
t = 0.55 eV, and a carrier density per unit cell n e ~ 10~ 3 
(which, according to Fig.^ puts EuBg far away from the 
WS) . If we incorporate such values in our previous re- 
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tal melts into a polaronic Fermi liquid or a metallic fer- 
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netic metal while at low temperatures, T < Tq, it per- 
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and the hopping energy t, and should apply to a broad 
class of materials that are described by the diluted DEM. 
In particular, our results provide a good description of 
the physics observed experimentally in EuBg. 
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